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Phase structure of a three-diensional Yukawa model 

Hidenori Sonoda 
Physics Department, Kobe University, Kobe 657-8501, Japan 

We use the method of the exact renormalization group to study the renormalization group 
flows of an O(N) invariant Yukawa model in three dimensional Euclidean space consisting of 
one real scalar and N real spinor fields. We obtain a phase structure similar to that of the N- 
vector model with cubic anisotropy, possessing a region of parameters exhibiting a first order 
transition. The particular case with one real fermion (N=l) belongs to the same universality 
class as the Wess-Zumino model with one supersymmetry. For the critical exponents of the 
Wilson-Fisher type fixed points, our 1-loop approximations are generally consistent with the 
results of previous studies. 

(N 

J>V §1. Introduction 

The purpose of this paper is to study a model of a real scalar field interacting 
\ with an arbitrary number of real spinor fields in three dimensional space-time. As 

is well known, the properties of a spinor such as the number of independent com- 
ponents depend highly on the dimensionality of space-timeP In three dimensional 
space-time, spinors transform under SL(2,R), and the real and imaginary parts do 
not mix. Under dimensional reduction, a Majorana spinor (or equivalently a Weyl 
spinor and its complex conjugate) in four dimensional space-time gives a complex 
spinor in three dimensions, which is a pair of real spinors. Hence, a single real spinor 
• is intrinsically an object in three dimensional space-time. Our model with an odd 

£> . number of spinors cannot be obtained from a four dimensional model by dimen- 

sional reduction. Therefore, we need a formalism that caters specifically to the three 
■ dimensional space-time. 

We study our model using the framework of the exact renormalization group 
£Sj | (ERG) upon which reviews from various viewpoints are available pW3 ERG is 

known for its successes in non-perturbative applications to many different areas of 
physics such as critical phenomena, gauge theories, lattice models, non-relativistic 
few-body problems, strongly correlated electrons, and quantum gravity.^ 

In this paper we are particularly interested in obtaining the phase structure of 
our model. To obtain the phase structure it is sufficient to derive the RG flows. To 
simplify this task, we restrict our model so that it is accessible from the Gaussian 
fixed point. With the O(N) symmetry among N real fermions, this choice leaves 
only three parameters to consider. Since we are mainly interested in understanding 
the general features of the model, we do not attempt to take advantage of the full 
potential of ERG. We solve the ERG differential equation for the Wilson action 
perturbatively at 1-loop level to derive the RG flows of the parameters. 

The primary finding of this paper is the phase structure similar to that of the N- 
vector model with cubic anisotropy.^ We find two non-trivial fixed points (besides 
the fixed point of the </> theory). One is the Wilson-Fisher type with only one 
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relevant parameter P3 For N = 1, this is the same as the fixed point of the M = 1 
supersymmetric Wess-Zumino model. For N even, this is the fixed point of the 
Gross-Neveu model with N/2 complex fermionsP^ The other fixed point has two 
relevant directions, and its existence implies a region of the parameter space that 
exhibits a first order transition. 

Our work has been partially motivated by a recent work on the M = 1 supersym- 
metric Wess-Zumino model in three dimensions™ where a non-trivial Wilson-Fisher 
type fixed point has been discovered. Our model includes the supersymmetric model 
as a subset, and we expect to obtain the same fixed point without imposing super- 
symmetry. If TV is even, the model can be rewritten for N/2 complex fermions, and 
this class of models has been studied with ERG in Ref . [T7|) , where the flow equation 
is solved numerically with the initial condition corresponding to the Gross-Neveu 
model. 

The paper is organized as follows. In £}2]we introduce our model and discuss its 
symmetry. In $3]we define the Wilson action and its parameters. In $5] we derive RG 
flows and find fixed points. In £j5]we discuss the nature of the phase transition shown 
by the model. In £j6] (for N = 1) and £J7] (for N even) we compare our results with 
those of previous studies. Finally, in $8] we conclude the paper with remarks. The 
extensive appendices provide technical details that make the paper self-contained. 

Throughout the paper we work in three dimensional Euclidean space, and we 
use the following notation for momentum integrals: 




§2. Yukawa model 

We consider a model whose classical lagrangian is given by 

c = \^ • + + \x j b ■ vx 1 + ^<4* V + ^ (2,1) 

where is a real scalar, and x 1 (I = lj • • • ,N) are real spinors. (See Appendix lAl 
for the corresponding lagrangian in Minkowski space.) We denote 

X = X T °y (2-2) 

We adopt the Einstein convention for summation over the repeated index /. 
The lagrangian is invariant under the following Z2 transformation 

</>(x) — > -<P(-x), x\x) — > ix\-x). (2-3) 

This invariance forbids the mass term xJx 1 ■ Depending on how this discrete sym- 
metry is realized, we expect two phases: 

1. Z2 exact — The expectation value of the scalar vanishes ((f)) = 0, and the 
fermions stay massless. 

2. Z2 spontaneously broken — ((f)) 7^ 0, and the fermions become massive. 
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In both phases we expect the 0(N) symmetry among N real fermions are unbroken. 

For N = 1, we expect that the model belongs to the same universality class as 
the N = 1 Wess-Zumino model whose classical lagrangian is given by 

£wz = \ ((V^) 2 + xa ■ V X ) + 9<f>lxX + 4- (0 2 - v 2 f (2-4) 



where 

m 2 = ~g 2 v 2 . (2-5) 

The Wess-Zumino model is a subset of the Yukawa model, satisfying the relation 

A = 3g 2 . (2-6) 

Hence, the critical exponents of the Wess-Zumino model must be the same as those 
of the Yukawa model. 

For N > 1, the Yukawa model also belongs to the same universality class as the 
three dimensional Gross-Neveu model given by 

£oK = I^.V*<-^(ixV) 2 . (2-7) 

We expect again that the critical exponents are the same as those of the Yukawa 
model. (See Fig.[TJ) For N even, the model contains N/2 complex fermions; defining 
complex fermions 

¥ = {x 2j ~ l + ix 2j ) , ff* = ^ (x*- 1 - (2-8) 
for j = 1, • • • , Yi we can write 

C GN = ¥a ■ VV j - ^ (W) 2 . (2-9) 
(See appendix [X] for complex fermions.) 

§3. Wilson action 

We construct a Wilson action Sa [<^ ; X 1 ] that depends on a UV cutoff A in such a 
way that its physics is A independent.— J The Wilson action is split into two parts: 

S A &, X 1 } = S F , A [<t>, X 1 } + S IA [<P, X 1 } ■ (3d) 
The free part is given by 

K {^4>{-p){p 2 + m 2 )(/)(p) + \pc\-V)W- ffX J (p)J • ( 3 '2) 



The cutoff function K(x) is positive, 1 at x 2 = 0, and decays fast enough for x 2 > 1. 
To make the physics independent of A, we impose that the interaction part Si a 
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m cr <g) N=l 



Fig. 1. Expected RG flows for the Yukawa theory with N real fermions. The Wess-Zumino (iV=l) 
and Gross-Neveu models [N > 1) belong to the same universality class as the Yukawa model. 



satisfy the ERG differential equation^ 

p)5(f>(p) } 



A JLq u i] = l [ A ^ A ^ f 5Sl ' A 6Sl > A i 6 Sl 
dA °lAi<P> XI 2 J p p 2 + m 2 \ 6<j){- P ) 5<P{p) 8^{-p)c 



(3-3) 



1 



TV (-i)p- a \ —j- rSi,A ■ Si,a , Tt \ + 7^77 \ S i,A~ 



2 J p p 2 [SX T (-P) ' ' Sx J (p) 5x J (-p) ' Sx J (p)j 

where 

A( q )^-2q 2 A K{q) . (3 . 4 ) 

To determine Si, a uniquely, we must introduce two additional conditions P^ 1 '^ 1 
1. UV renormalizability — We impose that the theory becomes the free mass- 
less theory at short distances. We demand the following asymptotic conditions: 

Si,a - J d 3 x (^,j/y~(V(£) 2 + mljy^cj) 2 + z^uv^a • Vy 7 

+_L<rtW~xV + ^uvt^) (3-5) 



2^ * uv 4! 

where the terms with higher powers of fields vanish, and all the coefficients are 

l uv 



constant except for mf JV which has a part linear in A and another linear in 



\nA//j,. 

2. Introduction of couplings A, g — We expand the Wilson action in powers 
of fields to obtain 

Si,a = j ^u 2 (A;p, -p) (f>(p)(f){-p) 
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Pl,P2,P3 



-^U 4 (A;p 1 ,p 2 ,P3,-pi ~ P2 ~ P3) 4>(Pl) ( t>(P2)(p(P3)4 ) (P4) 



+ 



p,k 



= G(A;p, k)4>(k)^(-p - k^ip) + 



(3-6) 



Choose a finite renormalization scale [i. We then impose 



u 2 (/x;0,0) 
" «a(A»;0,0) 



rip- 



9m 

9 n 2 (/i;p, -p) 



m 2 =0 
m 2 =0 



p 2 =m 2 =0 

*(/i;0) 



n 4 (^0,0,0,0) 

G(m;0,0) 



m 2 =0 
m 2 =0 
m 2 =0 



0, 
0, 
0, 
0. 



(3-7) 



-A. 
-0- 



The first four are normalization conditions, and the last two introduce coupling 
constants A &: g. Note that the squared mass parameter m 2 is introduced 
through 

The conditions (|3-4ti3-7j) determine Sa uniquely as a functional of (j) & x; Sa 
depends on the parameters m 2 ,X,g and the mass scales A, [i. Sa also depends on 
a particular choice of the cutoff function K, but it can be shown formally that the 
dependence can be absorbed by the normalization of parameters and fields. (For 
example, see Appendix B.2 of Ref. flQ|) .) 

The beta functions and anomalous dimensions are obtained as the \i dependence 
of the Wilson action JED MS 



d 



(3-8) 



The operators O m ,0\,O g are defined by 
8 



Ox 



dm 2 



Sa 



d_ 

d\ 
d_ 

dg 



K (p/A) (1-K (p/A)) 1 f SS A SS a | 5 2 S a 

{p 2 + m 2 ) 2 2 \ S(j>(p) 5<j){-p) 54>{p)5<p(-p) 

Sa , 



Sa- 



(3-9) 

(3-10) 
(3-11) 



These generate infinitesimal changes of the parameters m 2 , A, g, respectively. Denot- 
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ing the A independent correlation functions by brackets, we obtain 

f (O m 0( Pl )---x h (qi)---) m2Xg = -d m 2{<t){pi)---X h (qi)---) m 2 Xg , 

(OxHPi)---X h (Qi)---) m ,' Xg = -dx(4>(pi)---X h (qi)---) m . Xg , (3-12) 
{O g <P{pi)---X h {qi)---) m , Xg = -d g {4>{pi)---X h {qi)---) m , Xg , 

where the dots stand for a string of elementary fields <fr and x 1 ■ The operators A/^ 
and M x are the equation-of-motion operators defined by 

N^-jKWA^^tfM-VtfMj^ . (344) 

where 

Uln _ , Mj 1-K(p/A) 5S I)A ,„ n K > 

M W = 0(p) + p2+m2 ^(-p) ' (3 ' 15) 

[x 7 ] (p) = x 7 (p) + ^y^ H)?- °s¥Fp) Si ' a ■ (3 ' 16) 

They count the number of fields in the correlation functions: 
J <A^0(p 1 )---^ n )x /l ((Zi)---X /2fc fefc)) m2iAi5 = ™(0(Pi)---X 7l (<Zi)---) m2iAig , 

\ <AA x ^i)---^„)x /l ('7i)---x 72fc fe fc )> m2 ' A9 = 2/c<0( Pl )---x 7l (gi)---> m2iAi9 • 

(3-17) 

The meaning of (j3-8|) is clear. Its correlation with a product of elementary fields 
gives the RG equation: 

+ ^ + ^| + ^4) • • • Xh{qi) ■ ■ - 

= (Vfy + n xlx ) ((j)(pi) ■ ■ ■ x h {qi) ■ ■ -) m 2 iAi5 , (3-18) 

where n^,n x are the number of 0's and x' s m the correlator. 

In order to compute the beta functions at 1-loop, we need to compute Sj^a at 
1-loop, in particular the coefficients 112,114, z,G. The results are given in Appendix 
iBl Taking the \x derivatives, we obtain the following results: 

^ = -4{^ + 2I 6) + l(I 3 + 2I 5 )}, (3-20) 

/x " 2 /5 + ^iV /7 " ^3 (/3 + 2/5) ' (3 ' 21) 
7* = ^U3 + 2/5), (3-22) 
2 2 1 r 
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where the integrals I's are defined in terms of the cutoff function K and its derivative 
A in Appendix [El 

§4. RG equations and fixed points 

We have obtained 1-loop beta functions and anomalous dimensions. To obtain 
RG flows that describe the phase structure of the Yukawa model, we need to rescale 
both space and fields so that the renormalization scale [i is fixed under the RG flows. 
Due to this resclaing, the beta functions acquire contribution from the engineering 
dimensions of the parameters. Calling as m? , ^ as A, and as g, the flow 
equations for these dimensionless parameters become 



dt 
dg 2 
dt 



9 2 + 2<7# 



9 • 



(4-1) 



These equations are valid in a neighborhood of the origin m 2 = g 2 = A = 0, called 
the Gaussian fixed point. As has been explained in Jfl we restrict ourselves only to 
the region of parameters accessible from the Gaussian fixed point. In other words, 
we only follow the flows that originate from the origin. All the non-trivial fixed 
points we will discuss in this section are accessible from the Gaussian fixed point. 
At 1-loop, using the results of the previous section, we obtain 



f t m 2 

d \ 
dt A 



{2 - X\h ~ 9 2 \ (h + 2/ 5 )} m 2 + A±Ji - g 2 2I 2 , 
^(A + ( 5 Va) (\-\-(9 2 )) , 



(4-2) 



where the integrals J's are defined in terms of the cutoff function K in Appendix W\ 
gl and A/ by 



1 

j_ 

A/ 



= 1 (/ 3 + 2/5) + 2V (/ 5 + 2/ ( 



12 

3 T 
2^5 



and \±{g 2 ) by 



(4-3) 



\ + (g 2 ) + \.{g 2 ) 
\+(g 2 )\-(g 2 ) 



A / {l- 5 2 i(/ 3 + 2/ 5 )} 

N J 



-A/5 4 ^/ 7 . 



(4-4) 



With the convention A + (<? 2 ) > A_(<7 2 ), we obtain 



A ± (S 2 ) = ^ 



l-5 2 3(/ 3 + 2/ 5 ) 



(4-5) 



Similarly, the anomalous dimensions are given by 



70(5 2 ) 
lx(9 2 ) 



g 2 ±(I 3 + 2I 5 ) 



, . 2 

9 47v 



2 1 T 
— 15 



(4-6) 
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Fig. 2. Schematic RG flows. G for the Gaussian, I for the Ising, WF for the Wilson-Fisher, and 
B for the bicritical fixed points. The flow of m 2 is suppressed. We have a continuous phase 
transition in Region 1, and a first order transition in Region 2. 



which are independent of A. 

The RG flows are given schematically in Fig.[2j For any N = 1, 2, • • • , the flows 
have four fixed points. 

1. Gaussian — This exists by construction; in this paper we only study the RG 
flows out of this fixed point. At the Gaussian fixed point, all parameters vanish: 



m 



X = g = 0. 



(4-7) 



2. Ising — With g = 0, the fermions decouple, and we get a 4 theory. Its fixed 
point corresponds to the critical Ising model. 



rn 



A 

9 2 



A/ 
0. 



2-Aj 



h 

6/5-/4 



(4-8) 



The small deviation Am 2 = m? — mf has the scale dimension 

VE 



■2 -H —P n 



dm 2 



(4-9) 



Similarly, AX = A — A/ has the scale dimension —1. Am 2 and g 2 are the two 
relevant parameters at this fixed point. 
3. Wilson-Fisher — This is the most stable fixed point given by 



m 
A 

g 2 



m; 
A*- 

9*: 



-i-(A* + ^-^ 2 2/ 2 ) , 



VE+ 



(4-10) 



where A*+ is given explicitly by ()4-5p . The scale dimension of Am 2 is given by 

y E+ = 2-K + I -^-gl\{h + 2h) . (4-11) 
2 o 
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Table I. Notation for critical exponents 



parameter 


Gauss 


Ising 


Wilson-Fisher 


Bicritical 


Am 2 


2 


Vi > o 


VE+ > 


VE- > 


A\ 


1 


-1 


Vx+ < o 


vx- = -y\+ > o 




1 

2 


i 

2 


-1 


-l 



Am 2 is the only relevant parameter. A\ is irrelevant with the scale dimension 



1 

a7 



A*_) < 0. 



(4-12) 



4. Bicritical — This has two relevant parameters, and we call this a bicritical 
fixed point. 

m 2 _ = -^(A* + i- 5 2 2/ 2 ) , 



m 
A 



(4-13) 



5 2 =5*, 



where tje-, the scale dimension of Am 2 , is given by 



VE- 



(4-14) 



2 "*6 

Unlike the Wilson-Fisher fixed point, AX is relevant with the scale dimension 

>0. (4-15) 



y\- 



-y\- 



The critical exponents are universal, and accordingly they do not depend on 
what cutoff function K we use to formulate the ERG differential equations!^ (See 
also Appendix B.2 of Ref. I10p .) Approximate solutions of the ERG differential 
equations, whether they are perturbative or non-perturbative, lose universality, and 
the critical exponents become dependent on the choice of a cutoff function. The 
cutoff dependence is an inevitable artifact of the introduction of an approximation. 

To compute the critical exponents numerically, we have used a particular class 
of K (See Fig. [TO] in Appendix [E]): 



K{x) 



1 

l-x 2 





(x 2 < a 2 ) 
(a 2 < x 2 < 1) 
(x 2 > 1) . 



(4-16) 



The advantage of this cutoff function is that the necessary 1-loop integrals, listed 
in Appendix [El can be performed analytically. The cutoff function interpolates 
two popular choices, one at a = (Litim's cutoff 22 ') and another at a = 1 (sharp 
cutoff 23 '). The dependence of the critical exponents on the parameter a indicates 
the accurracy of the numerical values obtained. 

Let us plot the a-dependence of the critical exponents for N = 1. (Fig. We 
find that the exponents do not have much a-dependence for < a < 0.8. ue+ takes 
a maximum value 1.39 at a ~ 0.73, y\ + a maximum —0.71 at a = 0.98, and yE- a 
maximum 1.82 at a = 0. Since y\_ = —y\+, y\- takes a minimum 0.71 at a = 0.98. 
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0.2 0.3 0.4 0.5 0.0 0.7 



Fig. 3. Dependence of yE±,yx± on a for N = 1. We note j/a+ = —yx-- 

We take these extremum values as our numerical estimates for the exponents.^ 
Note that the following relations are satisfied: 



VE+ > y\+, HE- > V\- 



(4-17) 



Hence, at the bicritical fixed point, Am 2 is more relevant than AX. We next plot 
the a-dependence of the anomalous dimensions (Fig. 2]) : 



ri<p = 27^(5*), Vx = 2 7x(s*) • 



(4-18) 



At 1-loop, the Wilson-Fisher and bicritical fixed points share the same anomalous 
dimensions. Again, the dependence on a is mild in the region < a < 0.8. The 
extremum values at a = are 



V(p = Vx = 33 ~ ' 15 ' 



(4-19) 




Fig. 4. Dependence of J70,r? x on a for N = 1. Supersymmetry would imply % = rj x — 3 — 2y_B+. 
(See 31) 
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§5. Phase transitions 

Given A & g 2 , we examine the dependence of the model on the squared mass 
parameter m 2 in this section. 

Suppose (A,g 2 ) is in Region 1, surrounded by four RG flows G-I, I-WF, B-WF, 
and G-B in Fig. [2J Then the two-dimensional flow is always attracted to the Wilson- 
Fisher fixed point. Let m 2 r (X,g 2 ) be the value of m 2 such that the RG flow starting 
from (X,g 2 ,m 2 = m 2 r (\, g 2 )) flows to the Wilson-Fisher fixed point as t — > oo: 



< 



A 


t— »oo 






t— »oo 




fJ 2 




gl, 


2 


t—¥QO 


2 






mi 



+ ■ 



Since ue+ > 0, the deviation Am 2 = m 2 — m 2 + is relevant, and grows along the RG 
flow. Hence, the model exhibits a continuous phase transition at m 2 = m 2 r (X, g 2 ). 
We expect that the Z2 is exact for m 2 > m 2 r and broken for m 2 < m 2 r . 

Let us now suppose (X,g 2 ) is in Region 2, to the left of the RG flow G-B in 
Fig. [2j The two-dimensional flow has no fixed point to reach, and this implies the 
existence of a transition point m 2 r (\,g 2 ) that exhibits a first order phase transition. 
As for Region 1, we expect that the Z2 is exact for m 2 > m 2 r and broken for 
m 2 < m 2 r . In Appendix [C] we compute m 2 r and m 2 r by solving the 1-loop RG 
equations analytically. (For the first order transition, only the infinitesimal region 
below G-B is considered.) 

The phase structure of the Yukawa model discussed above is similar to that of 
the N-vector model with cubic anisotropy in 3 dimensions.^ (See also §5.8.5 of 
Ref. [25]) . for example.) The classical lagrangian of the model is given by 

N . N N 

£ = \Y.{ ( v ^) 2 + ™ 2 ^) + k E ti) 2 + 1 E > ( 5 - 2 ) 
1=1 1=1 ' 1=1 

where N > 2. The phase structure depends on N, but it is similar to the phase 
structure of the Yukawa model found here with the presence of a region of first order 
transitions. (Fig. [5]) 

§6. Comparison of the N = 1 case with the Wess-Zumino model 

The N = 1 case is particularly interesting since the model belongs to the same 
universality class as the N = 1 supersymmetric Wess-Zumino model. The Wess- 
Zumino model has recently been studied with ERG™ and we would like to make 
sure that our results are compatible. 

The Wilson-Fisher fixed point is characterized by three critical exponents: ue+, 
rjtf,, and r\ x . Supersymmetry implies 

V<f> = Vx- 

As Fig. H] shows, this relation is satisfied at the extremum a = 0, and reasonably 



12 



H. Sonoda 




N=2, 3 




N=5, 6, 



Fig. 5. RG flows of the N-vector model with cubic anisotropy in 3 dimensions — G, H, I stand 
for Gauss, Heisenberg, and Ising, respectively. First order transitions take place in the regions 
below G-I and on the left of G-C (N < 4) or G-H (JV > 4). H and C merge at N = 4. 

Table II. Numerical estimates for the critical exponents at the N = 1 Wilson-Fisher fixed point 





VE+ 


yx+ 


V4> 


Vx 


3 - 2y E+ 


Yukawa 


1.39 


-0.71 


0.15 


0.15 


0.22 


Wess-ZumincP 1 


1.406 


-0.756 


0.188 


0.188 


0.188 



satisfied for a G [0,0.8]. An interesting relation 

3 - 2y E+ = r) (j) = r) x (6-2) 

found in Ref. [16|) is also reasonably satisfied for a G [0.6, 0.8]. (Fig.HJ) Considering the 
crudeness of our 1-loop approximations, the agreement is satisfactory. In Appendix 
|D]we apply our perturbative ERG formalism directly to the WZ model and compute 
the critical exponents. The agreement with Ref. 116]) improves slightly. 

§7. Comparison with the Gross-Neveu model 

Let us consider the large N limit of the RG equations of sect. 01 which gives 

(7-1) 



VE+ 



1, 
1, 



VE- 

Vx 



3 ' 

0. 



The Gross-Neveu model with ^ complex fermions has been studied with ERG in 
[T7|) . where the RG flow of the Yukawa model was run numerically with the initial 
condition corresponding to the Gross-Neveu model. In TABLE IIIII the results of 
Ref. fT7|) are compared with ours for ue+, evaluated at a peak near a = 1, and 77^, 
evaluated at a = 0. (Fig. [6]) 

In Ref. [26]) . the critical exponents have been calculated to order 



jj as follows: 



VE+ 

V<t> 



16 J_ 

wr 

• 



(7-2) 



(We have replaced N, the number of complex fermions, in Ref. |26|) by where 
is the number of real fermions.) Our results reproduce the correct large N limit, but 
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Table III. Numerical estimates of critical exponents for large N. Our results are compared with 
those from Ref. 1 1 Tp . (In ll7|l . the notations v = l/yE+, r\ a — 77^ are used.) 





N 


= 4 


N 


= 6 


N 


= 8 


N = 


= 24 




VE+ 


V4> 


VE+ 


V4> 


VE+ 


V* 


VE+ 




ours 


1.16 


0.42 


1.11 


0.52 


1.08 


0.59 


1.03 


0.81 


Ref.O 


1.041 


0.561 


0.961 


0.710 


0.990 


0.789 


0.978 


0.936 



not the jr corrections. The most serious failure of our 1-loop approximation is the 
wrong sign for the correction to ue+] in comparison the correct sign has been 
obtained inll7p. 

We note that a first order transition was found for N = 2 in Ref. IT7|) . This can 
be explained if we assume that the transition point of this model belongs to Region 
1 of the (A,c/ 2 ) parameter space. (Fig. EJ) 




Fig. 7. The flow (left) starting from the transition point <?o,tr of the N = 2 Gross-Neveu model 
does not flow to a fixed point. This is compared with the flow (right) if the transition were 
continuous. 
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§8. Concluding remarks 

In this paper we have applied the ERG formalism perturbatively to the Yukawa 
model in three dimensions to obtain the RG flows. We have found a phase struc- 
ture similar to that of the N-vector model with cubic anisotropy. Accordingly, the 
spontaneous breaking of the Z2 symmetry of the Yukawa model can take place ei- 
ther at a continuous or at a first order transition point. The existence of a domain 
of parameters for the first order transition should be verified further, perhaps, by 
studying the effective potential with ERG. 

The Yukawa model with one real spinor includes the N = 1 Wess-Zumino model 
as a subset, and the fixed point of the latter is inevitably that of the former. At the 
Wilson-Fisher type fixed point of the Yukawa model, the only relevant parameter 
preserves supersymmetry; hence, the long distance behavior of an almost critical 
Yukawa model is the same as a massive Wess-Zumino model with its supersymmetry 
either exact or broken spontaneously. Emergence of supersymmetry at a critical 
point of a non-super symmetric theory has also been found for the Yukawa model 
with complex scalars and spinors P-^ 1 

The way we apply ERG perturbatively, we can only study the RG flows of UV 
renormalizable theories. This constraint still leaves a wide class of models as an 
object of study, and we advocate further perturbative applications of ERG for its 
simplicity and good cost-performance ratio. 
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Appendix A 

Spinors in D = 3 

For the reader's convenience, we summarize the salient features of spinors in 
three dimensional Minkowski and Euclidean spaces. 

A.l. Minkowski space 

We can choose the gamma matrices pure imaginary: 

7° = cr y , 7 1 =m a: , j 2 = ia z . (A-l) 
Under a Lorentz transformation, a two-component spinor x transforms as 

X^A X , (A-2) 

where A G SL(2,R) is a 2-by-2 real matrix with determinant 1. Since A is real, we 
can assume x to be real. We define 

X = X T o- y (A-3) 
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which transforms as 

X — > XA~ 1 , (A-4) 

since 

a y A T a y = A~ l . (A G SL(2, R)) (A-5) 
The lagrangian of the Yukawa model is given by 

C = \d^d^ - + - g^-x'x 1 ~ ^4> 4 • (A-6) 

This is invariant under the parity transformation defined by 



< 



4>(x,y,t) — > -4>{-x,y,t) , 

x(x,y,t) — > -i^x{-x,y,t) = a x x{-x,y,t), (A-7) 
, X(x,y,t) — > x{-x,y,t)i^ = x(-x,y,t)(-a x ) . 



The mass term \X is forbidden by this invariance. 

A Dirac spinor is a linear combination of two real spinors: 



# = -^ty + #'). (A-8) 



We find 

a, = a/t^/ = 

We obtain 

# {i^j^d^i -M)V = - {^{i^d^ - M)ij> + ^'(i^d^ - M)y/) . (A- 10) 

A Dirac spinor can be obtained by the dimensional reduction of a Weyl (or Majorana) 
spinor in (3 + l)-dim Minkowski space. 

A. 2. Euclidean space 
We can choose 

7 = a (A-ll) 

corresponding to the spin \ representation, familiar from non-relativistic quantum 
mechanics. Under a space rotation, a two-component spinor transforms as 

tfj^Uip, (A- 12) 

where U G 577(2). tp = tp T a y transforms as 

$ — »• V^" 1 . (A-13) 

Hence, the lagrangian 

C = ^(i>v ■ V-0 + rwiptpj (A- 14) 

is invariant under SU{2). 
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We call ip a real fermion, and call a linear combination of two real fermions ip, tp' 

1 



(A-15) 



a complex fermion. Denoting \P = -^=(ip — itp'), we can write the lagrangian for W as 



C = V + m)>f. 

•Z' and <F are two independent 2-component Grassmann fields. 



(A-16) 



Defining 



Appendix B 

Wilson action at 1-loop 



d 



(B-l) 



the 1-loop vertices are calculated as follows: 
1. scalar 2-point 



u 2 (A;p, -p) = ^ 



A(q) 



1-K(q/A) 1 /• 1 - K(q) \ 
q2(q2 _|_ m 2^ n I n A 



+9' 



pf_ 
2 



(q+pf 



2 N 



n 



l-K(q/A)l-K((q+p)/A) | 1 f( l-K(q) } 



q 2 J 



+ / 1 {(1 - K (q/A)) (1-K ((q - p)/A)) - (1 - K (q/A)) 2 } 
+£V-2(/l-/x) / 1^)(1-K( 9 )) , 



where the constant c 2 is given by 

1 f 1 

6 



C2 



(l-K(g))(Z\(g)-4(g)). 



(B-2) 



(B-3) 



2. fermion 2-point 



z 2 (/l;p 2 



iV 



(q+pf 



q 2 + m 2 



+ 



q 2 + m 2 



2 



i-^(g) V 

<7 2 ) 



l-K((q+p)/A)l-K(q/A) _ f 1 - K (q/A) y \ qp + h 



J ( P 



, (B-4) 



(q+p) 2 

where the constant k 2 is given by 

k 2 = \J^(l- K(q)) {-2(1 - K(q)) + A(q)} . (B-5) 
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3. Yukawa coupling 



G(A;0,0) = -g + ^ 



{1-K{q/A)f 1 [(l-K(q)f 
q 2 (q 2 + m 2 ) fi 



(B-6) 



4. scalar 4-point 



~ 6 N /^{a-^^M)) 4 -^ 1 "^)) 4 } • (B-7) 



The asymptotic behaviors as 7l — )• oo are given by 



1 

12 



w 2 (^l;M,-M) 



l 



G(i;0,0) 



i« 4 (A-o,o,o,o) 



A J_ 

2l2 



zA(g) _ m_ 2 /" 1 - 
<? 2 



2 A-p, f A(q)(l-K(q)) 



A 2 



^g 2 -2 

H p 

15 2 



1 



+ C 2 



(l-K(g)) 2 



+ fe 2 



1 g 3 f(l-K(q)f 



g 

VI n ^yfA, 

A_3 ^_ [(l-K(q)) 2 
A 2fiA 



+ 



f(l-K(q)Y 

N jiA 



(B-8) 
(B-9) 
(B40) 



(B41) 



Appendix C 

Analytic calculations ofm 2 r (X,g 2 ) andm 2 r (X,g 2 ) 



The 1-loop RG equations can be solved analytically. Hence, given arbitrary 
(X,g 2 ) in Region 1, we can compute the critical value m 2 r . Similarly, given (\,g 2 ) 
in Region 2 near the line GB, we can compute the transition value m 2 r . In this 
appendix, we first solve the flow equations, and then calculate m 2 r and m 2 r . 

C.l. Solving the flow equations 

The flow equations derived in 2] are given by 



ig 2 (t) = 5 2 (i)(l-4f 



dt' 



A(t) = j-{\ + {g 2 (t))-\{t)}{\{t)-\_{g 2 {t))} , 



(C-l) 
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X G 



X X- 



Fig. 8. RG flows 



where X±(g 2 ) (A + > A_) are defined by 



\ + (g 2 )\-(g 2 j - 

The solution for g 2 {t) is given by 
x{t) =g 2 {t)e- 1 

We choose < g 2 (0) < g 2 so that 



X + (g 2 ) + X-(g 2 ) = X I {l-g 2 l(I 3 + 2I 5 )} 

■ 6 T 

TV 7 ?- 



1 + 



1 e* 



x(t) ^°°0. 



Defining 



we obtain 



A(t) 



A, 



/(*) 

— = k ■ x(t) y *- , 



m - % 

where A; is a constant determined by the initial condition, and 

A*_i- — A* 



y\- = 



Aj 



> 0. 



In Fig. [HJ the three regions Rf, R l , R2 are defined by the behavior of / 



J 9* 



R 



+ . A, 
1 
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k is positive in Rf, but negative in R± and R%- In R2, the region of a first order 
transition, we find / — > —00 (hence A — > —00) at a finite t. 
The flow equation for the squared mass is given by 

jm 2 {t) = {2 + Ag 2 {t) + BX(t)} m 2 (t) + Cg 2 (t) + DX(t) , (C-9) 



where 



A = --(I 3 + 2I 5 ) 



B = ~\h, C = -2I 2 , D = ~h. 



(C40) 



This is solved as 



2t- [ dt' (Ag 2 {t') + BX(t')) 
Jo 

+ f dt' {Cg 2 (t') + DX(t')) exp -2t' - f dt" (Ag 2 (t") + BX(t")) 
Jo Jo 



m 2 (t) = m 2 (0) 
■t' 



.(Oil) 



C.2. Critical squared mass m 2 r (X,g 2 ) 

In region i?i = R^ U -Rj - , the critical value of m 2 (0) = m 2 is obtained as a 
function of g 2 (0) = g 2 and A(0) = A by the condition 



'" -> mi. 



m 2 (t) 



(C-12) 



This implies that the left-hand side of ()C-lip vanishes in the limit t —> +00. Hence, 
we obtain 



(ft(C# 2 (i) + £>A(i))exp -2i- / dt' (Ag 2 (t') + BX{t')) 

L 



(C-13) 



Using the analytic solution for g 2 (i) and A(t), this can be rewritten as 



m 2 r (A,5 2 ) 



2 f 1 dx C + Df(x) 



1 + 4 f 1 -! 



exp 



9 2 ['^(A + Bfix')) 

J x x 



(C-14) 

where /(x) is the same / as before except it is regarded now as a function of x = 
g 2 (t)/g 2 , and it is explicitly given by 



A _ A* 
<F 9: 



^ + ( 4? - A ) ■ % 



A _ A*- 



I A»+ A 



(C-15) 



C.3. First order transition point m 2 r (X,g 2 ) 

We consider a region directly below the trajectory GB, which is a line given by 

" 2 X 



9- =- _ 

5* ^* 



(C-16) 
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For g < g 2 and a positive infinitesimal e, 

X = ^( 9 2 + e)<^g 2 (C47) 

gives a point in R2 just below the trajectory GB. The squared mass at the first order 
transition is obtained as 



m 2 r ( ~2~ {9 2 +e) ,5 2 J - ( hrg 2 ,g 2 J + ^-e -^-m 2 cr (X, g 2 ) 



9* / \ 9* / 9* 



9; 



(C-18) 



where 



m cr ( —9 ,9 = VE-m*- — 2 / eta 

9* y M-2! 



= Tpi: /„ 2 d2; » ( C49 ) 

and the derivative is given by 



1 



R „2 

a 2 / \ 2^ 

-ml. (A. ) 



— m 2 cr (X,g 2 ) 



B 9* 2(2 2 

V .., ,,, ~ 72 m cr ( 9 —,9 

9 



Appendix D 

ERG for the M = 1 Wess-Zumino model in D = 3 

To study the three dimensional M = 1 Wess-Zumino model, it is the most 
convenient if we preserve supersymmetry manifestly by linearizing it with the help 
of a real auxiliary field F. The corresponding classical lagrangian is 

Cwz = ^V</> • V<t> + l -xd -V X + \F 2 + g jiF~ (0 2 - v 2 ) + <f^xx\ • (D-l) 

Integrating over F, we reproduce the classical action given by (|2-4[) in sect. [2j The 
imaginary i is necessary for the Euclidean space. The classical lagrangian is invariant 
under the following N = 1 supersymmetry transformation: 

H = xf, 

= (*-V0-iF)f, (D-2) 
<5iF = Vx-^, 
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where £ is an arbitrary constant spinor. Besides super symmetry, the theory is in- 
variant under the Z2 transformation: 



( F(x) 



-<f>(-x) , 
ix(-x) , 
F(-x). 



(D-3) 



We construct a Wilson action that is invariant under the linearized supersymmetry 
and Z2 transformation. (The Z2 invariance forbids the supersymmetric mass term 
iF<f> + ±xx-) 

The Wilson action is split into the free and interaction parts: 

Sa = ~\J p K (p/A) ( p2 ^- p ^ + K-ti B • + F (~P) F ip)) + S ^ > ( D4 ) 

where Sj yi satisfies the ERG differential equation 



1 



A(p/A) 



A M Sl ' A ~2j p 



-Tr(-i)p- a 



+ 



S 2 Sr A 



54>(-p) 6<t>(p) ' 5<j){-p)54>{p) 

7 „ „ V 1 



^x(-p) 



SV/i • <SV. 



+ 



5/, 



5 



+ 



$x(p) $x{-p) ' $x(p) 



SF(-p)5F(p) 5F(-p)5F(p) J _ 
For UV renormalizability, we impose the following asymptotic behavior: 



(D-5) 



A— >oo 



-J, 



1 



^v^((V0) 2 + x^-Vx + F 2 ) 



(D-6) 



where and gt/y are constants, and (<7f 2 )r/v is linear in A (Note that in the 
common language the coefficient of iF(0) is linearly divergent; supersymmetry is 
consistent with this divergence.) The parameters g, v 2 are introduced through the 
expansion: 

Si,a = ~ J z^A-p) 1 - (p 2 <P(-p)Hp) + X(-P)i° ■ PX(p) + F(-p)F(p)) 



J ■ 

Jp,k 



G(A;p, q) - {iF(-p - q)4>{p)4>{q) + (f>(~p - q)x{p)x{l)) 



+(GV 2 )(A)- i F(0), 



where we impose 



* 2 (a*;0) = 0, 

G(/x;0,0) = g, 
{ (GV 2 )(v) = gv\ 



(D-7) 



(D-8) 
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at an arbitrary renormalization scale \i. 

The dependence of Sa on \x is given by the RG equation 

- h^Sa = -pd g S A - p v 2d v 2S A + iM , (D-9) 

where N is the equation-of-motion operator that counts the number of fields. The 
1-loop calculations give the coefficients as 



fa = -fih + ^v 2 (I 6 + \h) , (D-10) 



I 4u J 5i 



where the integrals are defined in Appendix [El Unlike the N = 1 Wess-Zumino 
model in 4 dimensions or its dimensionally reduced N = 2 model in 3 dimensions, 
there is no non-renormalization theorem for this model. 

Redefining g 2 / \i by g 2 and v 2 / ' \i by v 2 , the RG flows are given by 



fg 2 = g 2 + 2g/3, 
f t v 2 = v 2 + V 



! 2 ;7"' (D-ll) 



Besides the Gaussian fixed point g 2 = v 2 = 0, the flows have a non-trivial fixed point 
at 



2 _ 1 1 
1 



J 2 ~ / 



(D-12) 



The anomalous dimension of the elementary fields is given by 



r, = 2 1 {gl) = - ^ , (D-13) 
4 h + jh 



and the scale dimension of v 2 — v 2 is 



1 h + \h 



1 h + ih 



These 1-loop results confirm the sum rule found in Ref. [R 

2y E + -n = Z (D-15) 

which is valid for any choice of the cutoff function in our case. Using the particular 
cutoff function given in Appendix [Ej we obtain the plot in Fig. [9) The exponent 
He — 1-43 hardly depends on the cutoff parameter, and it agrees a little better with 
Ue+ = 1.408 from Ref. [TFj) than ue+ — 1.38 of the N = 1 Yukawa model. 
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Fig. 9. The exponent yE — 1.43 for the Wess-Zumino model is compared with yE+ — 1-38 of the 
Yukawa model. 



Appendix E 

Integrals of a cutoff function 



We define the following integrals: 



h = . 


lq g 2 ' 


h = 




r ^(<?) 2 


h = 


h = > 


r a-K(q)f 

lq q 4 ' 


h = 


h = . 


r (i-K(q)) 4 

»9 q 4 ■ 





r A(q)(l-K(q)) 
Jq q 1 

Jq Tj 4 ' 



For a particular choice of the cutoff function (see Fig. [10 



K(q;a) 



1 

1-a 2 



for < q 2 < a 2 

for a 2 < q 2 < 1 
2 



for K. 
the above integrals can be evaluated easily as follows: 



h(a) 
Ha) 
h(a) 
Ha) 



1 
1 
1 

2^ 
1 

2^ 



2 l+q+q 2 

3 1+a ' 

4 1+a+a 2 

3 (l+a) 2 (l-a) ' 

4 l+2a 

3 (1+a) 2 ' 

8 5+20a+29a 2 +16a 3 



Ha) 
Ha) 
Ha) 



1 

27T 3 
1 

27T 3 
1 

27T 3 



2 2a 3 +4a 2 +6a+3 
15 (1+a) 2 : 
2 

1+a' 

6+18a+16a 2 
5(l+a) 3 ' 



(E-l) 



(E-2) 



(E-3) 
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